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We study the extrapolation ol nuclear shell structure to the region of superheavy nuclei in self- 
consistent mean-field models - the Skyrme-Hartree-Fock approach and the relativistic mean-field 
model - using a large number of parameterizations which give similar resuhs for stable nuclei 
but differ in detail. Results obtained with the Folded- Yukawa potential which is widely used in 
macroscopic-macroscopic models are shown for comparison. We focus on differences in the isospin de- 
pendence of the spin-orbit interaction and the effective mass between the models and their infiuence 
on single-particle spectra. The predictive power of the mean-field models concerning single-particle 
spectra is discussed for the examples of ^"^Pb and the spin-orbit splittings of selected neutron and 
proton levels in ^®0, ^'^^Sn and ^"^Pb. While all relativistic models give a reasonable description 
of spin-orbit splittings, all Skyrme interactions show a wrong trend with mass number. The spin- 
orbit splitting of heavy nuclei might be overestimated by 40% — 80%, which exposes a fundamental 
deficiency of the current non-relativistic models. In most cases the occurrence of spherical shell 
closures is found to be nucleon-number dependent. Spherical doubly-magic superheavy nuclei are 
found at fg|ll4, ^72120, or fg2l26 depending on the parameterization. The Z = 114 proton sheU 
closure, which is related to a large spin-orbit splitting of proton 2/ states, is predicted only by forces 
which by far overestimate the proton spin-orbit splitting in ^''^Pb. The Z — 120 and N — 172 shell 
closures predicted by the relativistic models and some Skyrme interactions are found to be related 
to a central depression of the nuclear density distribution. This effect cannot appear in macroscopic- 
microscopic models or semi-classical approaches like the extended Thomas-Fermi-Strutinski integral 
approach which have a limited freedom for the density distribution only. In summary, our findings 
give a strong argument for 172 120 to be the next spherical doubly-magic superheavy nucleus. 
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I. INTRODUCTION 

The extrapolation of nuclear shell structure to super- 
heavy systems has been discussed since the early days of 
the shell correction method ||l|-^, when spherical proton 
shell closures at Z — 114 and Z = 126 and a spherical 
neutron shell closure at TV = 184 were predicted. Shell 
effects are crucial for the stability of superheavy nuclei 
which by definition have a negligible liquid-drop fission 
barrier. Recent experimental progress allowed the syn- 
thesis of three new superheavy elements with Z = 110- 
112 [p[-pO|, but these nuclides are believed to be well 
deformed. The experimental data on these nuclei and 
their decay products - a-decay half-lives and Qa val- 
ues - agree with the theoretical prediction ]ll|-[l6[| of a 
deformed neutron shell at A^ = 162 which has a signif- 
icant stabilizing effect [ |lO|jl7| ]. The experimental proof 
of the deformed shell by a measurement of the deforma- 
tion is beyond the current experimental possibilities. As 
a first step in this direction the ground-state deforma- 
tion of ^^'*Noio2 was deduced from its ground-state ro- 
tational band in a recent experiment |1^ . The ultimate 
goal is to reach the expected island of spherical doubly- 
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magic superheavy nuclei. More refined parameterizations 
of macroscopic-microscopic models [13-16 
older finding that it is located around jgl 114 
nuclei, although even heavier than the heaviest nuclides 
known so far, are expected to have much longer half-lives 
due to the stabilizing effect of the spherical shell closure 
which significantly increases the fission barriers |l9-p2]. 

Although modern macroscopic-microscopic models 
quite successfully describe the bulk properties of known 
nuclei throughout the chart of nuclei, their parameteri- 
zation needs preconceived knowledge about the density 
distribution and the nuclear potentials which fades away 
when going to the limits of stability. Like the mean- 
field models based on the shell correction method, self- 
consistent mean-field models have been used for the in- 
vestigation of superheavy nuclei from the earliest param- 
eterizations |^,Q to the most recent ones psj-p^. 

In two previous articles we have discussed the occur- 
rence of spherical |^ and deformed ||3^ shell closures 
in superheavy nuclei for a large number of parameteriza- 
tions of self-consistent nuclear structure models, namely 
the Skyrme-Hartree-Fock (SHF) approac h [p3 , and the 
relativistic mean-field (RMF) model |35| - |37[ . Spherical 



proton shell closures are predicted for Z — llA, Z = 120 
and Z = 126, depending on the parameterization, while 
neutron shell closures occur at A^ = 172 and N — 184 re- 
spectively. Only one parameterization - the Skyrme in- 
teraction SkI4 - confirms the prediction of macroscopic- 
microscopic models for a doubly magic ^11114, other pa- 
rameterizations - the Skyrme forces SkM* and SkP - 
predict 5g4l26, while yet others - the Skyrme interac- 
tion SkI3 and most of the relativistic forces - give a new 
alternative with ^72120. Several interactions predict no 
doubly magic spherical superheavy nucleus at all. In self- 
consistent models, the proton and neutron shells strongly 
affect each other |Q . Small details of the shell structure 
have a strong influence on the potential energy surfaces 
of superheavy nuclei in the vicinity of the ground-state 
deformation, leading to dramatic differences in the fis- 
sion barrier heights and therefore in the fission half-lives, 
while the predictions of different models and forces are 
similar at large deformations [p3[ . 

Superheavy nuclei differ from stable nuclei by their 
larger charge and mass numbers. The strong Coulomb 
potential induces significant changes in the proton shell 
structure: single-particle states with large angular mo- 
mentum and small overlap with the nuclear center only 
are lowered compared to small-j states, see Figs. 1-2 of 
Ref. JS^I and the discussion therein. While this effect oc- 
curs already in non-sclf-consistent models, polarization 
effects of the density distribution due to the high charge 
number can be described in self-consistent models only. 
The Coulomb interaction pushes protons to larger radii, 
which changes the density distribution and the single- 
particle potentials of both protons and neutrons in a 
complicated manner. On the other hand, the large mass 
number of superheavy nuclei leads to a high average den- 
sity of single-particle levels. Therefore the search for shell 
effects in superheavy nuclei probes the detailed relations 
among the single-particle states with extremely high sen- 
sitivity. 

The question arises which features of the effective 
mean-field models are most decisive for the single-particle 
structure. The three most crucial ingredients in this re- 
spect are: first the effective nucleon mass and its radial 
dependence which determines the level density near the 
Fermi surface, second the spin-orbit potential which de- 
termines the energetic distance of the spin-orbit partners, 
and third the density dependence of potential and effec- 
tive mass which has an influence on the relative position 
of the states. We perform here a comparison of various 
parameterizations from SHF as well as RMF with empha- 
sis on their spin-orbit properties. The effective masses 
(with one exception) are comparable in all forces. The 
density dependences are similar amongst the SHF forces 
and amongst the RMF forces, but differ signiflcantly be- 
tween SHF and RMF. The largest variations in the sam- 
ple of parameterizations occurs indeed for the spin-orbit 
part of the forces where we have three classes, the stan- 
dard SHF models, SHF with extended spin-orbit forces 
(SkI3 and SkI4), and the RMF models. The present pa- 



per concentrates predominantly on this given variation 
of the spin-orbit force. It is the aim of this paper to ex- 
plain the contradicting results of self-consistent models 
mentioned above and to find the most reliable prediction 
for the next spherical doubly-magic superheavy nucleus. 
In Sect. H the properties of the mean-field models and 
the parameterizations used is discussed. In Sect. HI the 



details of the spin-orbit interaction and the differences 
between the various models used are explained, while 
Sect. 1^ discusses briefly the relation between effective 
mass and average density of single-particle levels. In 
Sect. M we compare the predictions of the various mean- 
field models with known single-particle energies in ^"^Pb 
and experimental spin-orbit splittings in ^^O, ^^^Sn and 
^"^Pb and study the shell structure of the potential spher- 
ical doubly-magic nuclei ^g4ll4, ^y^l^O, and 184I26 and 
the predicted nucleon- number dependence of the Z — 120 
proton shell and the N = 172 neutron shell in some de- 
tail. Sect. VI summarizes our findings. In an appendix 
we present the details of the mean-field and pairing mod- 
els necessary for our discussion. 



II. THE FRAMEWORK 

The Skyrme force was originally designed as an ef- 
fective two-body interaction for self-consistent nuclear 
structure calculations. It has the technical advantage 
that the exchange terms in the Hartree-Fock equations 
have the same form as the direct terms and therefore the 
numerical solution of the Skyrme-Hartree-Fock equations 
is as simple as in case of the Hartree approach, while 
the solution of the Hartree-Fock equations using finite- 
range forces like the Gogny force [g8| is a numerically 
challenging task. The total binding energy can be for- 
mulated in terms of an energy functional which depends 
on local densities and currents only, see Appendix H. 
This links the Skyrme-Hartree-Fock model to the effec- 
tive energy functional theory in the Kohn-Sham approach 
which was originally developed for many-electron sys- 
tems. The Hohenberg-Kohn theorem ^^ states that the 
non-degenerate ground-state energy of a many-Fermion 
system with local two-body interactions is a unique func- 
tional of the local density only. The Kohn-Sham scheme 
[Q relies on the Hohenberg-Kohn theorem but keeps the 
full dependence on the single-particle wavefunctions for 
the kinetic energy which allows to preserve the full shell 
structure while employing for the rest rather simple func- 
tional in local-density approximation. This point of view 
can be carried over to the case of nuclei where, however, 
the non-local two-body interaction requires an extension 
of the energy functional by a dependence on other densi- 
ties and currents, e.g., the spin-orbit current. In any case, 
there is no need for a fundamental two-body force in an 
effective many-body theory, but one can start from an 
effective energy functional which is formulated directly 
at the level of one-body densities and currents (see, e.g., 



| pl[ and references therein). 

The relativistic mean-field model can be seen from the 
same point of view as a relativistic generalization of the 
non-relativistic models using a finite-range interaction 
formulated in terms of effective mcsonic fields. Relativis- 
tic kinematics plays no role in nuclear structure physics, 
but the RMF naturally describes the spin-orbit interac- 
tion in nuclei, which is a relativistic effect that has to 
be added phenomenologically in n on-relativistic models. 
This will be discussed in Sect. 



Ill in more detail. 



For both SHE and RMF there are numerous parame- 
terizations in the literature. We select here a few typical 
samples of comparable (high) quality, mostly from recent 
fits. For the non-relativistic Skyrme-Hartree-Fock calcu- 
lations we consider the Skyrme forces SkM* Ml, SkP 
H, SLy6, SLy7 |4|g|, Skll, SkI3, and SkI4 Q. For 
the RMF we consider NL3 0, NL-Z, Q, and two com- 
pletely new forces, NL-Z2 and NL-VTl. All forces are 
developed through fits to given nuclear data, but with 
different bias. Of course, the basic ground-state proper- 
ties of spherical nuclei (energy, radius) are always well 
reproduced. Small variations appear with respect to fur- 
ther demands. The parameterization SkM* is the oldest 
in the list here. It was the first Skyrme force with accept- 
able incompressibility as well as fission properties and re- 
mains up to date a reliable parameterization in several 
respects. The Skyrme force SkP was developed around 
the same time with the aim to allow the simultaneous 
description of mean-field and pairing channel. Moreover, 
it was decided here to use effective mass m* /m = 1.0. 
(Mind that all other forces in our sample have smaller 
effective masses around 0.6 < m* /m < 0.8). The forces 
SLy6 and SLy7 stem from a series of fits where it was suc- 
cessfully attempted to cover properties of neutron mat- 
ter together with normal nuclear ground-state proper- 
ties. In SLy6 the contribution of the kinetic terms of 
the Skyrme force to the spin-orbit potential is discarded, 
which is common practice for nearly all Skyrme parame- 
terizations, e.g. SkM* and the Sklx forces in the sample 
here. SLy7 is fitted exactly in the same way as SLy6, but 
these additional contributions to the spin-orbit force are 
considered, see the discussion in Sect. Ill A for details. 
The forces Skll, SkI3 and SkM stem from a recent series 
of fits along the strategy of E^] where additionally key 
features of the nuclear charge formfactor were included 
providing information on the nuclear surface thickness. 
For these, furthermore, information from exotic nuclei 
was taken into account in order to better determine the 
isotopic parameters. The force Skll is a fit within the 
standard parameterization of the Skyrme forces. This 
performs very well in all respects, except for the isotopic 
trends of the charge radii in the lead region. To cover 
these data, one needs to extend the spin-orbit functional 
by complementing it with an additional isovector degree 
of freedom Efl] as will be discussed in Sect. [II A in more 
detail. SkI3 uses a fixed isovector part built in analogy to 
the RMF, whereas SkM was fitted allowing free variation 
of the isovector spin-orbit force. The modified spin-orbit 



force has a strong effect on the spectral distribution in 
heavy nuclei and thus even more influence for the predic- 
tions of shell closures in the region of superheavy nuclei. 
The forces headed by "NL" belong to the domain of 
the RMF model. The parameterizations NL-Z, NL-Z2, 
and NL3 use the standard nonlinear ansatz for the RMF 
model, whereas NL-VTl additionally considers a tensor 
coupling of the vector mesons. The parameterization NL- 
Z H^ aims at a best fit to nuclear ground-state proper- 
ties along the strategy of |49| . It is a re-fit of the popular 
force NLl with a microscopic treatment of the correction 
for spurious center-of-mass motion. NL-Z2 and NL-VTl 
are new parameterizations developed for the purpose of 
these studies to match exactly the same enlarged set of 
data including information on exotic nuclei like the Sklx 
Skyrme forces. This should allow better comparison be- 
tween the RMF and the Skyrme model. The force NL3, 
finally, results from a recent fit including neutron rms 
radii. It gives a good description of both nuclear ground 
states and giant resonances. Details of the RMF La- 
grangian and the actual parameterizations are discussed 
in Appendix A 2 . 



The nuclear matter properties of the forces are sum- 
marized in Table ffl. These are to be considered mainly as 
extrapolations from finite nuclei to the infinite system. 
There a few exceptions because in some cases the one or 
the other nuclear matter property has entered as a con- 
straint into the fit. These cases are: the effective mass 
m*/m = 1 for SkP, the compressibility Koo = 230 MeV 
and asymmetry coefficient agym = 32.0 MeV for the SLyx 
forces, and the sum-rule enhancement factor k = 0.25 in 
case of the SLyx and Skia; forces. Table || shows that 
most Skyrme forces share the basic nuclear matter prop- 
erties close to the phenomenological values like binding 
energy per nucleon E/A « — 16 MeV, equilibrium den- 
sity po « 0.16 fm"^, incompressibility K = 210±30MeV 
M 1 asymmetry energy 30 MeV < agym < 32 MeV, and a 
low sum-rule enhancement factor < k < 0.25. A phe- 
nomenological value for the effective mass of m* /m w 0.8 
can be drawn from the position of the giant quadrupole 
resonance in heavy nuclei |5l| . And we see that the mean 
field results for the effective mass vary in a wide range 
0.58 < m* /m < 1.0 about this value. This is a bit disqui- 
eting because the effective mass is a feature which has a 
strong impact on spectral properties, influencing, in turn, 
the predictions for superheavy nuclei. 

The nuclear matter properties of the relativistic pa- 
rameterizations differ significantly from those of Skyrme 
forces. E/A is usually slightly larger and po somewhat 
smaller than the values for Skyrme interactions. The pre- 
dictions for the incompressibility K differ systematically 
from those of the nonrelativistic models, in case of NL3 
it is somewhat larger, in case of the other RMF forces 
smaller than the average result for Skyrme forces. But 
all parametrizations stay within the accepted bounds of 
this rather uncertain quantity. The asymmetry coeffi- 
cient and the sum-rule enhancement factor are substan- 
tially larger than in case of the Skyrme forces. But all 



TABLE I. Compilation of nuclear matter properties for the parameter 
equilibrium energy per nucleon and density, Koo the compression modulus, 
(note that we provide two values for the relativistic models where the value 
the other at fc = fl 



sets used in this study. E/A and po denote the 
m* /m the effective mass in units of the free mass 
in brackets is m* /mikp) at the Fermi surface and 
flaym the asymmetry coefficient, and k, the sum-rule enhancement factor. 



Force 


E/A [MeV] 


po [fm-^] 


K^ [MeV] 


m* jm 


tlsym 


K 


SkP 


-16.04 


0.163 


202 


1.000 


30.0 


0.35 


SkM* 


-16.01 


0.160 


217 


0.789 


30.0 


0.53 


SLy6 


-15.92 


0.159 


230 


0.690 


32.0 


0.25 


SLy7 


-15.90 


0.158 


230 


0.688 


32.0 


0.25 


Skll 


-15.93 


0.160 


243 


0.693 


37.5 


0.25 


SkI3 


-15.96 


0.158 


258 


0.577 


34.8 


0.25 


SkI4 


-15.92 


0.160 


248 


0.650 


29.5 


0.25 


NL3 


-16.24 


0.148 


272 


0.595 (0.659) 


37.4 


0.68 (0.53) 


NL-Z 


-16.18 


0.151 


173 


0.583 (0.648) 


41.7 


0.72 (0.55) 


NL-Z2 


-16.07 


0.151 


172 


0.583 (0.648) 


39.0 


0.72 (0.55) 


NL-VTl 


-16.10 


0.150 


179 


0.600 (0.663) 


39.0 


0.66 (0.51) 



RMF forces agree in their rather low value for the effec- 
tive mass 0.58 < m* /m < 0.6. It is to be noted, however, 
that the effective mass in RMF depends on the momen- 
tum as 



m*{kF) I fm*{0)y , fkp'^ ^ 




m*(0) 



0.08 



(1) 



where m*{0) is the value at A: = usually handled as 
effective mass in the RMF and where we assumed in the 
second step a typical kp ~ 1.35/fm. Table | thus shows 
two values for m* /m in case of the RMF, at momentum 
zero and in brackets the more relevant value m* /m(kp) 
at the Fermi surface. The latter value is larger by about 
10% and comes visibly closer to the results for the Skyrme 
forces. 

In view of the application to superheavy nuclei, it is 
worthwhile to check the performance of all these forces 
in our sample with respect to already known superheavy 
nuclei. This was done in Ref. ||3^. It turns out, that SkI3, 
SkM, and the relativistic forces perform best in that re- 
spect, although it is to be mentioned that all relativistic 
forces show a wrong isotopic trend, see ||3^ for details. 
It is noteworthy that the extended Skyrme functionals 
SkI3 and SkM perform much better in the region of su- 
perheavy nuclei than the Skyrme parameterizations with 
the standard spin-orbit interaction. This indicates that 
an extended spin-orbit interaction is an essential ingre- 
dient for the description of heavy systems. 

In both SHF and RMF the pairing correlations are 
treated in the BC S scheme using a delta pairing force, 
see Appendix A3 for details. 

The numerical procedure solves the coupled SHF and 
RMF equations on a grid in coordinate space with the 
damped gradient iteration method |53] . The codes for the 
solution of both SHF and RMF models have been imple- 



mented in a common programming environment sharing 
all the crucial basic routines. 



III. SPIN-ORBIT INTERACTION IN NUCLEAR 
MEAN-FIELD MODELS 

A. The spin-orbit field 

The spin-orbit interaction is an essential ingredient 
of every model dealing with nuclear shell structure 
to explain the shell closures of heavy nuclei beyond 
N ^ Z ^20 [|3|J54). It was aheady noted in the first 
explorations with the modified oscillator model that dif- 
ferent fits of the spin-orbit coupling constant lead to con- 
tradicting predictions for the next major shell closures in 
superheavy nuclei |55| . 

The spin-orbit interaction emerges naturally in rela- 
tivistic models and the explanation of the large spin- 
orbit splitting in nuclei was one of the first prominent 
successes of the relativistic mean-field approach |5^] . The 
spin-orbit potential can be deduced in the non-relativistic 
limit of the RMF and is given up to order v^/c^ by Pq] 



"3 



he 



{2m q + Sq 



Ty^'<^« 



Va 



(2) 



where S and V are the scalar and vector potentials re- 
spectively, see appendix A 2 for details. While the usual 



potential is given by the sum of the large negative scalar 
potential S and the large positive vector potential V 
which cancel nearly to give the usual shell-model poten- 
tial, the difi^erence of scalar and vector potential enters 
the expression for the spin-orbit field, explaining its large 
strength. The occurrence of the derivative of the fields 
in (0) indicates that the spin-orbit field is peaked in the 
nuclear surface region and that its strength will depend 
on the surface thickness of the particular nucleus. 



To compare with the corresponding expression for 
Skyrme interactions, one has to evaluate (^ in local den- 
sity approximation 



\\A(™f) 
9 



he 



(2m, ~Cp- Cp^y 



{CVp + C'Vpq 



(3) 



where C = C^ 



a 



Cp and C" = 2Cp are combina- 
tions of RMF parameters with Ci — gf/rn^. The isospin 
dependency of the spin-orbit potential is rather weak for 
typical RMF parameterizations which give C" « 0.1 x C. 
In the framework of non-relativistic models the zero- 
range two-body spin-orbit interaction proposed by Bell 
and Skyrme p7|j5q | is widely used. Examples are all stan- 
dard Skyrme interactions like SkM*, SkP, the SLyx forces 
or Skll and other non-relativistic effective interactions 
like the Gogny force p8|. The corresponding spin-orbit 
potential Wg is given by 



Wfd)=64(Vp + Vp,) 



(4) 



There are two fundamental differences between the rela- 
tivistic and non-relativistic expressions for the spin-orbit 
potential: the isospin dependence and the missing den- 
sity dependence in case of the non-relativistic models. 

When deriving the single-particle Hamiltonian from an 
underlying Skyrme force there appears an additional con- 
tribution to the spin-orbit field which arises from the 
momentum-dependent terms in the two-body Skyrme 
force 



W(J) = &4(Vp + Vp,) + ciJ - c'lJ, 



(5) 



The calculation of the spin-orbit current J is somewhat 
cumbersome in deformed codes and its contribution to 
the total binding energy rather small. Therefore the J- 
dependent terms in (0) are discarded in most parameteri- 
zations of the Skyrme interaction and ^ is used instead. 
SkP and SLy7 are two exceptions in this investigation. 

In the Hohenberg-Kohn-Sham interpretation of the 
Skyrme interaction outlined above, there is no need for 
an underlying two-body force, but one can start from an 
effective energy functional which is formulated directly at 
the level of local one-body densities and currents. This 
relaxes the fixed isotopic mix (]^ in the spin-orbit func- 
tional and allows more freedom for its parameterization 
which was used to complement the spin-orbit interaction 
by an explicit isovector degree of freedom in the fit of the 
extended Skyrme functionals SkI3 and SkI4 



9 



64 \/p + &4 \/pq 



(6) 



The additional isospin dcgrce-of-frecdom enables the re- 
production of the kink in the isotope shifts of charge 
mean-square radii in lead, which is not possible with stan- 
dard Skyrme forces employing (||) P6| , ^j60| , while the 
experimental data are reproduced by most RMF forces. 
The parameters 64 and 64 in SkI3 and SkI4 are adjusted 
to reproduce the spin-orbit splittings of protons and neu- 
trons in ^^O and the isotope shifts of charge mean-square 



radii in lead. As a result of the fit the approximate re- 
lation 64 » —6 4 emerges for SkM, see also Table H in 



Appendix A 1 



This means that for SkM the spin-orbit 
potential of one kind of nucleons depends mainly on the 
density profile of the other kind of nucleons. The force 
SkI3 was adjusted with the same fit strategy but with 
a fixed isovector part 64 = analogous to the RMF in 
the sense that the spin-orbit potentials of protons and 
neutrons are approximately equal. However, there re- 
main differences between SklB and the RMF: all RMF 
potentials have a finite range and the spin-orbit interac- 
tion has a small but non-zero isospin dependence and a 
strong density dependence. 



B. Spin— orbit splitting 

In non-relativistic models the spin-orbit term in the 
equation-of-motion of the radial wavefunctions in case of 
spherical symmetry is given by 



Wg,^[]k{jk + l)-ik{h + l) 



(t>k{r) 



(7) 



where Wq^r is the radial component of the spin-orbit 
potential and (j){r) the radial part of the single-particle 
wavefunction V'(r). For well-bound single-particle states, 
the radial wavefunctions 4)t±i/2 entering Eq. (M) are 
only slightly different. Therefore the contributions from 
the potential and the kinetic term can be neglected in 
very good approximation when calculating the spin-orbit 
splitting Aels — £^'+1/2 ~ ££-1/2 of two states (j)j with the 
same radial quantum number and orbital angular mo- 
mentum i but different j = ^ ± 1/2 



AeLS«47r / dr r T^,,,(£ + i) {l+ l)\(j,,^^j^\ 



i(^+h)\(^i-i/2? 



(8) 



The spin-orbit splitting scales with 2£ + 1 and depends 
sensitively on the overlap of the single-particle wavefunc- 
tions with Wq^r/T- The shape of Wq,rlr - which is usu- 
ally peaked at the nuclear surface - depends itself on the 
variation of the actual density distribution in the nucleus 
which changes going along isotopic or isotonic chains, es- 
pecially when the density distribution becomes diffuse 
going towards the drip-lines or when it develops a cen- 
tral depression - as happens in some superheavy nuclei, 
see Sect. VD|. 

Equation (^ holds as well for the non-self-consistent 
single-particle models which are used in the framework 
of macroscopic-microscopic models. There the spin-orbit 
potential W is assumed to be proportional to the gra- 
dient of the single-particle potential U . In the simplest 
case of the modified oscillator model - which was used 
in the first studies of the shell structure of superheavy 
nuclei pj|] - the spin-orbit potential W/r has no radial 



dependence, the amplitude of the spin-orbit spUtting is 
simply proportional to 2£ + 1, see |5^] for a detailed dis- 
cussion. In more refined single-particle models like the 
Folded- Yukawa model (FY) |6^ or Woods-Saxon model 
p2[ the spin-orbit potential is peaked at the nuclear sur- 
face like in the self-consistent models, see appendix A4 
for details. 



IV. EFFECTIVE MASS AND AVERAGE LEVEL 
DENSITY 

The average density of single-particle levels g{e) in the 
vicinity of the Fermi energy can be estimated using the 
Fermi gas model in a finite potential well. In case of 
non-relativistic particles one obtains |63] 



S2g{N,) = E{N, - 2) - E{N,) 



(10) 



fff' 



2m* 






(9) 



The relativistic generalization of formula (^ is simply 
obtained by inserting the effective mass m*{kp) at the 
Fermi surface, see eq. (|l|) and the values in brackets in 
table 1). 

The average level density rises linearly with particle 
number - the single-particle spectra of superheavy nuclei 
are therefore much denser than those of lighter stable nu- 
clei. This makes the shell structure of superheavy nuclei 
very sensitive to details of the spin-orbit interaction, dif- 
ferences of a few 100 keV in the spin-orbit splitting of 
two given orbitals can create or destroy shell closures. 

The level density depends linearly on the effective mass 
m* as well. This causes a dramatic difference when com- 
paring the predictions of interactions with small effective 
mass, e.g. SkI3 with m* /m = 0.574, and parameteriza- 
tions with large effective mass like SkP with m* /m — 1.0 
in the region of superheavy nuclei. As said before, a 
phenomenological value of m* /m k, 0.8 for the isoscalar 
effective mass can be determined from the position of the 
isoscalar quadrupole giant resonances which is just in be- 
tween the extremes spanned by our choice of mean field 
models. But a word of caution is in place here. The value 
of 0.8 is appropriate for the effective mass in the nuclear 
volume. But the value may be larger at the surface, or 
Fermi surface respectively p4| . This is, admittedly, a 
feature which is not yet built into nowadays mean field 
models. A thorough exploration of this aspect is a task 
for future reasearch. 



V. SPHERICAL MAGIC NUMBERS 

A. Relation of single-particle spectra and bulk 
properties 

At closed shells, one observes a sudden jump in the 
two-nucleon separation energies 5*2^ 



Nq and the number of the other kind of nucleons are 
assumed to be even. The two-nucleon separation energy 
is a better tool to quantify shell effects than the single- 
nucleon separation energy due to the absence of odd-even 
effects. It is a very good approximation for twice the 
negative Fermi energy 



S2q{N,Z)^-2\{N,Z) 



(11) 



In doubly-magic nuclei - in which the BCS pairing model 
breaks down - the Fermi energy is simply given by the 
single-particle energy of the last occupied state. Devi- 
ations between the calculated and experimental values 
for the single-particle energy of the last occupied state in 
doubly magic nuclei are therefore connected by (O) with 
an error in the two-nucleon separation energies below the 
shell closure. Although slightly influenced by pairing cor- 
relations, this holds in good approximation also for the 
first unoccupied state above the Fermi surface and the 
two-nucleon separation beyond the shell closure. 

The size of the gap in the single-particle spectrum is 
given by half the difference in Fermi energy when going 
from a closed shell nucleus to a nucleus with two addi- 
tional like nucleons. But from Eq. ( pi] ) it follows that 
this is in very good approximation equal to the shell gap 
52q, the second difference of the binding energy 



52q{Nq) = E{Nq + 2) - 2E{Nq) + E{Nq - 2) 
^-2[\{Nq + 2)-\{Nq)] , 



(12) 



which was used in |3l[| to quantify the magicity of a nu- 
cleus. Going away from closed shells, there is a non- 
negligible contribution from the residual pairing interac- 
tion, therefore 5*2^ and ^2g loose their direct relation to 
the single-particle levels. The two-nucleon gaps 52q rep- 
resent the size of the gap in the single-particle spectra, 
but they do not contain information about the actual 
location of the single-particle energies. 

Only interactions which reproduce the experimental 
values of the first single-particle state below and above 
the Fermi surface will give the correct binding energies 
around closed shell nuclei. This can be read the other 
way round as well: Only interactions which reproduce 
the binding energies around shell closures give a good 
description of at least the first single-particle state be- 
low and above the shell closure, but the bulk properties 
give no information on single-particle states away from 
the Fermi energy. This demonstrates nicely, however, 
that total binding energy and properties of single-particle 
states are connected in self-consistent mean-field models. 
This is very different in macroscopic-microscopic models 
where the bulk properties and single-particle spectra are 
described in separate models. 

One has to be careful when comparing experimen- 
tal and calculated single-particle spectra. Experimental 
single-particle energies of even-even nuclei are deduced 



from excitation energy measurements of adjacent odd- 
mass nuclei. The binding energy of odd-mass nuclei is 
affected by polarization effects induced by the odd nu- 
cleon, see |6^] for a discussion of these effects in the 
framework of the RMF. The polarization effects are im- 
portant for the comparison of calculated and experimen- 
tal single-particle energies. But they do not affect the 
relation between the single-particle spectra and the bulk 
properties in even-even nuclei discussed here. 



B. The single- particle spectra in known nuclei 

Before extrapolating the models to the regime of su- 
perheavy nuclei we want to test the predictive power 
of the mean- field models looking at ^°*Pb, the heav- 
iest known spherical doubly-magic nucleus. Figure |l| 
shows the single-particle spectra of ^°®Pb as obtained 
from spherical calculations with the mean-field forces as 
indicated. The upper panel shows the spectrum of the 
protons, the lower panel that of the neutrons. The exper- 
imental excitation energies in the neighboring odd nuclei 
are shown on the left side for comparison, the data are 
taken from JG^]. The gaps in the single-particle spectra 
at Z = 82 and N = 126 are clearly visible, but the forces 
obviously do not agree for this stable nucleus, which was 
used in the fit of all parameter s ets em ployed here. 

As already discussed in Sect. VA , the difference be- 
tween the calculated and experimental energies of the 
first single-particle state above and below the shell clo- 
sure refiects the quality of the description of the total 
binding energies in the vicinity of a shell closure. There 
are large differences between the forces in their predic- 
tions for states further away from the Fermi surface. The 
spectrum predicted by SkP is much too dense and the 
ordering of proton states below the Fermi surface not re- 
produced. A natural explanation for this might the too 
large effective mass of SkP, but one has to be careful: 
The effective mass determines the average level density 
only but not the level density in an actual nucleus. The 
difference in energies between the 2gg/2+ and lzii/2+ neu- 
tron states is, for example, by far too large when calcu- 
lated with SkP and SkM*, leading to a sub-shell closure 
at iV = 136 in contradiction to experimental data. In 
the RMF and extended Skyrme forces this difference is 
by far too small, NL3 predicts even a wrong ordering 
of these two levels. The relativistic forces and the rel- 
ativistic corrected Skyrme force SkIS overestimate the 
gap between the proton I/19/2- and 2/7/2- states above 
the Fermi surface which leads to a pronounced sub-shell 
closure at Z = 92 which again is in contradiction with 
experiment. 

The RMF models and the modern Skyrme forces with 
small effective mass push the IJ15/2- with an experi- 
mental single-particle energy of e = —2.51 MeV too much 
up in the spectrum, e.g. to e = -0.418 MeV in NL-VTl, 
while Skyrme forces with a large effective mass like SkM* 
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FIG. 1. Single-particle spectrum of the protons (upper 
panel) and neutrons (lower panel) in '^"^Pb calculated with 
the mean-field forces as indicated. 



and SkP work slightly better within this respect. The dif- 
ferences in average level density due to the actual value 
of the effective mass scale only the deviation from the ex- 
perimental value. States with large orbital angular mo- 
mentum systematically lie too high in the single-particle 
spectrum for all forces, see also the proton lii3/2+ state. 
As this problem appears for all parameterizations of both 
SHF and RMF models and for all nuclei throughout the 
chart of nuclei |£^,^, we conclude that this is not a 
problem of actual fits but it indicates the need for im- 
proved effective interactions beyond the current energy 
functionals. 

All forces have problems to reproduce the neutron 
single-particle energies below the Fermi energy as well. 
All relativistic forces and SkIS give a wrong level or- 
dering, the 2/5/2- state lies too low in energy in all 
cases. Standard Skyrme forces work slightly better in 
that respect, e.g. SkP predicts 2/5/2- to be the second- 
to-last state below the Fermi surface, but interchanges 
the 3p3/2- and 1^13/2+ states instead, the latter one is 
again pushed up too much in energy like all other states 
with large angular momentum. It is remarkable that the 
non-self-consistent FY model is the only one which re- 
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FIG. 2. Relative error {5eis - Sel^'^^)/5el'^^' in percent of 
the spin-orbit splitting of proton (upper panel) and neutron 
(lower panel) single-particle states close to the Fermi surface 
in ^^O, ^^^Sn, and ^°*Pb calculated from the mere mean-field 
single-particle energies with the parameterizations as indi- 
cated. 



produces the level ordering of all states in the vicinity of 
the Fermi energy for both protons and neutrons. Like 
the self-consistent models, however, it is not able to re- 
produce the values of the single-particle energies or even 
their relative distance. 

To conclude our findings so far: the comparison be- 
tween predictions of various current mean-field mod- 
els and experimental data shows that the models are 
not able to reproduce all details of experimental single- 
particle spectra and show additionally significant differ- 
ences among each other which are related to effective 
mass and details of the spin-orbit interaction. 

Shell closures of heavy nuclei are related to the spin- 
orbit splitting of states with large orbital angular mo- 
mentum. Therefore it is interesting to compare the pre- 



dictions of the mean-field models with experimental data 
on spin-orbit splittings in known nuclei. Figure |^ shows 
the relative errors in % of the spin-orbit splittings of neu- 
tron levels (lower panel) and proton levels (upper panel) 
near the Fermi surface in ^^O, ^■^^Sn and ^"^Pb. Nega- 
tive errors denote theoretical values which are too small. 
The spin-orbit splittings are calculated from the single- 
particle energies as they come out from a spherical mean- 
field calculation. As already mentioned, the experimental 
single-particle energies are measured as separation ener- 
gies between adjacent nuclei, where polarization effects 
have a visible influence. The error bars in Fig. rep- 
resent the uncertainty of the spin-orbit splittings due to 
polarization effects as they are found in |^ . 

All RMF forces reproduce the experimental spin-orbit 
splittings fairly well, although there are deviations up to 
20% which are scattered around zero. The errors from all 
RMF forces are similar and therefore it is likely that these 
errors represent the standard RMF Lagrangian, not spe- 
cific parameterizations. Although the tensor couplings 
of the vector mesons in NL-VTl change the relative dis- 
tance of the single-particle energies compared to NL-Z2, 
see Fig. y, they have no visible influence on the spin-orbit 
splittings compared to the standard Lagrangian. It is in- 
teresting that the errors of the spin-orbit splittings of the 
neutron ip and 2/ states in 208p|^ have the largest val- 
ues but different sign while ^^O and ^■^^Sn are described 
very well. There is only one splitting known for protons 
in ^"^Pb (if one excludes splittings across the Fermi sur- 
face which have a large theoretical uncertainty, see |]65| ) , 
so one has no information how the error depends on the 
angular momentum of the state as in the case of neu- 
trons. But, however, the RMF gives a very good overall 
description of spin-orbit splittings throughout the chart 
of nuclei without any free parameters adjusted to single- 
particle data. 

The reproduction of the experimental data with the 
Skyrme functionals is by far not as good as for the rela- 
tivistic models. There is a clear trend which is the same 
for all standard Skyrme forces: for neutrons the error 
of the Ip splitting in ^®0 has the smallest value, then 
comes the splitting of the 3p state in ^°*Pb, the 2d state 
in ^^^Sn and then the splitting of the 2/ state in ^°^Pb. 
Like in the case of the RMF, the splittings of the 2/ and 
3p neutron states in ^"^Pb are not reproduced with the 
same quality, the error for the 2/ state is always much 
larger compared to the experimental value than for the 
3p state. 

It is very unlucky that the parameters of the spin-orbit 
interaction in non-relativistic models are usually adjusted 
to data in ^^O, which is at the lower end of a systematic 
trend increasing with mass number. Choosing one or 
several heavier nucleus for the fit, however, does not cure 
the problem of the wrong trend, but it gives a better 
overall description of spin-orbit splittings as can be seen 
from SkP, which gives the best possible compromise for 
a standard Skyrme force: the differences between the 
data points are similar to those from the other standard 
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FIG. 3. Two-proton gap in the chain oi Z — 114 isotones 
calculated with the forces as indicated. 



Skyrine forces, but they are centered around zero. The 
other standard Skyrme forces, SkM*, SLy6, SLy7, and 
Skll, give similar predictions, with large errors for the 
2d states in ^^^Sn and the neutron 2/ and proton 2d 
state in ^''^Pb. 

The predictions of the extended Skyrme forces SkI3 
and SkM deviate significantly from both the standard 
Skyrme forces and the RMF. SkI3 gives bad results for 
neutrons and protons and shows surprising large differ- 
ences to the relativistic forces. This is somewhat un- 
expected because SkI3 was constructed with the isospin 
dependence of the spin-orbit force which appears in the 
relativistic models. This indicates that the isospin de- 
pendence is not the only important difference between 
the relativistic and non-relativistic models, density de- 
pendence or finite range of the RMF potentials might 
play a much larger role for the single-particle spectra. 
SkM, gives the best results for the neutrons of all non- 
relativistic models, but at the same time it gives also 
the worst description for the proton spin-orbit splittings 
among all interactions investigated here, the errors have 
values up to 80% for the 2d level in ^ospb. The predic- 
tions for heavy nuclei might be too large by a factor of 
nearly two, which makes the unique prediction of this 
force of a proton shell closure an Z = 114, caused by 
large spin-orbit splitting, not ve ry re liable. This will be 
discussed in more detail in Sect. VC. 



The Folded- Yukawa model shows a similar behavior as 
the SHF forces, but like in case of SkP the errors are 
scattered around zero. 



C. The shell structure of llflU 

The nucleus ill 114 is the "traditional" prediction for 
the spherical doubly magic superheavy nucleus P,p|j55| 
from macroscopic-microscopic models which was con- 
firmed in more recent models of this type |l4|-[l6[. As 
shown in |^,Q, most modern parameterizations of 
self-consistent models shift this property to larger pro- 
ton numbers and/or smaller neutron numbers, depend- 
ing on the parameterization. Only for the extended 




FIG. 4. Single-particle spectra of igflW for protons (top) 
and neutrons (bottom) at spherical shape for the mean-field 
forces as indicated. 



Skyrme functional SkM 1I4II4 remains the doubly spher- 
ical magic nucleus in the superheavy region. 

Figure y shows the two-proton shell gap S2p, the indi- 
cator for shell closures derived from total binding ener- 
gies, for the chain of Z = llA isotopes calculated with the 
mean-field forces as indicated. Only SkM predicts a shell 
closure for Z = 114, all other forces give rather small 62p- 
In contrast to the proton shell closures at higher charge 
numbers Z which will be discussed in the following, the 
Z = 114 shell is stable for varying neutron number. 

We want to see now how the different predictions for 
the shell gap 52p in the potential doubly-magic nucleus 
^^^114 are reflected in its single-particle spectra, see Fig- 
ure Q. The possible shell closure at Z = 114 is located 
between two spin-orbit coupled states, the 2/7/2- and 
2/5/2- levels. Additionally, the 1^13/2+ state which has 
a similar energy as the 2/ states has to be pushed down. 
Therefore it is immediately clear that Z — 114 is only 
magic in the case of a large amplitude of the spin-orbit 
splitting. A strong Z = 114 shell appears only for SkM, 
the force with the largest proton spin-orbit splitting in 
this nucleus of all forces under investigation. But it is 
to be remembered that SkM overestimates the spin-orbit 



splitting of the protons in ^°*Pb by 80%. This makes the 
prediction of a large spin-orbit splitting in ^^*114, leading 
to a strong shell closure, very doubtful. 

SkP, the force with effective mass m* /m = 1.0 and 
therefore a large density of single-particle levels shows 
no significant shell structure at the Fermi surface of the 
protons at all. For all other forces there is at least a sub- 
shell closure dX Z = 114. But only for SkI4 the gap in the 
single-proton spectrum is large enough to be interpreted 
as a major shell closure. For all standard Skyrme forces 
the l«i3/2+ state is located between the 2/ states, which 
significantly reduces the Z = 114 gap. 

In some of the other forces with smaller spin-orbit split- 
ting, like SkI3 and the RMF parameterizations, there is 
a gap in the spectrum at Z = 120 indicating the major 
shell closure of these forces, while in all Skyrme forces 
there appears a gap at Z — 126, hinting at another po- 
tential spherical magic proton number. But as we will 
see in what follows the gap at Z — 126 becomes smaller 
with increasing proton number and has disappeared for 
most of the forces when reaching this proton number. In 
the single-particle spectrum for the neutrons in ig4ll4 
the differences between the various mean-field forces are 
much smaller than for the protons. All forces show a 
gap in the single- neutron spectrum at A^ = 184, but for 
the relativistic parameterizations the amplitude of this 
gap is smaller than for the Skyrme forces and even de- 
creases with increasing effective mass. Therefore, in NL3 
(the RMF force with the largest effective mass) the major 
shell closure at iV = 184 has vanished. 

The single- particle spectra of both protons and neu- 
trons from the non-self-consistent FY model look very 
different compared to all self-consistent models. In par- 
ticular, the spin-orbit splitting of all proton states is 
much larger compared to all self-consistent models with 
the exception of SkI4. At the Fermi surface, the lii3/2+ 
proton state which is the last filled state in all standard 
Skyrme forces, is pushed down below the 2/7/2- state by 
the large spin-orbit splitting. This creates the large gap 
in the single-particle spectrum at Z = 114. 

Although the non-self-consistent FY model predicts 
N — 184 to be magic as well, the ordering of the neutron 
states below the A^ — 184 shell closure is very different. 
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FIG. 6. Single-proton levels in the vicinity of the 

Fermi energy for the isotopes oi Z — 120 (lower panel) and 
two-proton shell gap S2p (upper panel) versus the neutron 
number, computed with SkI3. The dotted line in the upper 
panel is twice the difference between the 3pi m- and 2/5/2- 
levels, the dashed line twice the difference between the lin /2+ 
and 2/g/2- levels. 



The large spin-orbit splitting in the FY model pushes 
the 3(i3/2+ state above the 4si/2+ state and the 8^5/2+ 
below the 2^^7/2+ state. Another difference to the self- 
consistent models is the large level density above the gap 
at A^ = 184. Three states with large angular momentum, 
i.e. 2/111/2- , 1/13/2- and 1^17/2+ are close together which 
explains that the maximum of the corresponding shell 
correction is shifted to nuclei with the somewhat smaller 
(and non-magic!) neutron number around A^ = 178 p6[. 
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FIG. 5. Two-proton gap in the chain oi Z — 120 isotones 
calculated with the parameterizations as indicated. 



D. The Z = 120 shell 

In self-consistent models, the occurrence of a spherical 
proton shell closure with given Z can change with varying 
neutron number A^, and similarly the neutron shell clo- 
sures can vary with changing proton numbers. While for 
light nuclei this happens only at the limits of stability, 
e.g. the vanishing of the A^ = 28 shell for proton num- 
bers Z < 20 which is hinted experimentally [ [70| [73| and 
predicted by self-consistent mean-field models 17 4y75[] . In 
the region of superheavy nuclei the nucleon number de- 
pendence of shell closures is a common feature in the 
predictions of self-consistent models fplUsT 
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FIG. 7. The same as in Fig. 0, but computed with SkI4. FIG. 8. The same as in Fig. M, but computed with SkP. 



The most important example is the spherical Z — 120 
shell, see Fig. Is] which shows the two-proton shell gap S2p 
of the Z = 120 isotones for some of the forces under in- 
vestigation. All parametrizations except SkM* and SkP 
predict a peak in the 52p a.t N = 172 which is followed 
by a steep decrease of S2p when going towards larger 
neutron numbers. The i52p are largest in the relativis- 
tic parametrizations and the extended Skyrme functional 
SkI3 with the RMF-like spin-orbit interaction, but even 
most of the standard Skyrme forces, i.e. those with small 
effective mass, show an enhanced S2p around N = 172 as 
well. 

To understand the origin of the neutron-number de- 
pendence. Fig. H shows the single-proton spectra (lower 
panel) and the corresponding S2p (upper panel) of the 
Z = 120 isotones calculated with SklS. The quantity of 
interest is the gap in the spectrum at Z ~ 120. First of 
all it is to be noted that the single-particle spectrum is 
indeed relatively dense. Therefore already minimal rel- 
ative changes of the proton levels produce a regime of 
higher level density at the proton Fermi surface around 
N — 184, the neutron number where the proton shell gap 
is lowest. The relative changes of the levels are due to 
changes in the amplitude of the spin-orbit splitting. The 
shell closure a.t Z = 120 can appear only when the spin- 
orbit splitting between the 2/ proton states below the 
Fermi energy and the 3p states above the Fermi energy 
is small. In nuclei for which the spin-orbit splitting of 
these levels is large, e.g. around N — 184, the gap in the 
single-particle spectrum ai Z — 120 vanishes. 



To demonstrate the relation between the shell gap cal- 
culated from total binding energies and the actual gap 
in the single-particle spectrum, in the upper panel of 
Fig. the difference in energy Ae between the last single- 
particle state below and the first state above the Fermi 
energy is shown with a dotted line. As can be clearly 
seen, Ae is always larger than i52p, showing that the shell 
gaps 52p calculated from total binding energies are influ- 
enced by the pairing, which smears out the shell effects. 

For SkM the spin-orbit splitting of the single-proton 
levels in superheavy nuclei is in general larger than for 
Skl3, see Fig. 0. Therefore the magic number Z = 114 
appears, corresponding to a large gap between the 2/ 
single-proton levels. Like for SkI3, the spin-orbit split- 
ting of the levels in the vicinity of the Fermi energy is 
largest around N — 184. While this effect weakens the 
shell gap at Z = 120 in Skl3 and SkM, it amplifies the 
gap in the single-proton spectrum at Z = 114 in SkM. 
The magic Z = 120 appears for SkM only for isotopes 
with relatively small spin-orbit splitting in the vicinity 
of the Fermi energy, i.e. at large neutron numbers. 

The single-particle spectra of the protons look very dif- 
ferent for forces with large effective mass, e.g. SkP, see 
Fig. 0. Owing to the large average level density at the 
Fermi surface there are no distinct shell effects at all for 
the Z — 120 isotopes. Additionally, there are only slight 
changes of the level structure with varying neutron num- 
ber N. This confirms our previous finding that a large 
effective mass washes out most of the shell structure in 
superheavy nuclei. In this case, the proton shell gap 62p 
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FIG. 9. The same as in Fig. H, but for ?ral20 
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FIG. 10. Density distribution (upper panel) and radial 
component of the spin-orbit potential (lower panel) of pro- 
tons (right), neutrons (left) for 17I12O, calculated with the 
forces a indicated. The total density is plotted in the up- 
per panels as well. The density distributions calculated from 
the single-particle wavefunctions as they come out in the FY 
model is drawn for comparison. All models except SkP show 
a central depression in the density distribution, which has a 
visible impact on the spin-orbit potential. 



and the last single-particle level below the Fermi energy 
and the first level above are in good agreement. 



E. The Shell Structure of ??il20 

The occurrence of the proton shell closure aX Z = 120 is 
coupled to at least a subshell closure at iV = 172. There- 
fore it is interesting to take a detailed look into the single- 
particle spectra of f 7I 120, which are shown in Fig. g. The 
upper panel shows the proton levels, the lower one shows 
the neutron levels. As already discussed in Sect. VD, 
the occurence of the shell closure a.t Z = 120 depends on 
the amplitude of the spin-orbit splitting of the 3p states 
above the Fermi level and the 2/ levels below the Fermi 
energy. It appears only when the level density at the 
Fermi energy is small and the spin-orbit splitting is weak, 
but this is the case for all forces under investigation ex- 
cept SkP and SkM*, the forces with the largest effective 
mass and therefore largest (average) level density. It has 
to be noted that for almost all forces this nucleus is lo- 
cated near the two-proton drip line since the first unoc- 
cupied proton level has a positive single-particle energy. 



The level ordering of the proton states above the Fermi 
level for the RMF forces NL-Z, NL-Z2 and NL-VTl is 
quite unusual, the 3p state with small total angular mo- 
mentum is located above the state with large angular 
momentum. This phenomenon is related to the unusual 
shape of the density distribution of this nucleus, see the 
upper panel of Fig. |l^. The large dip at the nuclear 
center, where the density is reduced to 2/3 of its nuclear 
matter value, leads to a region around r « 3 f m where the 
spin-orbit potential has the opposite sign, see the lower 
panel of Fig. nO. Therefore, for j states with large oc- 
cupation probability in this region the amplitude of the 
spin-orbit splitting is dramatically reduced or even has 
the opposite sign as it is the case for NL-Z, NL-Z2 and 
NL-VTl. Additionally, this density distribution strongly 
affects the shape of the single-particle potentials, which 
are reduced at the nuclear center by approximately the 
same factor as the density. Orbitals with large angular 
momentum, e.g. the li states, are pushed down in the 
spectrum compared to states with rather small angular 
momentum like the 3p states. This leads to a completely 
different level ordering above the Z — 120 proton shell in 
case of the RMF forces. 

The same effect occurs in the neutron spectrum as well. 
The level ordering of the 3d states is reversed for the 
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calculated with Skl3. The probability distribution is shown for the state with larger total angular momentum only. 



RMF forces, see the lower panel of Fig. y. Again, for 
SkP, the force which gives the less pronounced dip of 
the density distribution, the spin-orbit splitting of the 
3d neutron states is largest. States with large angular 
momentum and therefore small overlap with the center 
of the nucleus, i.e. the 2g or Ij states, show the common 
spin-orbit splitting. 

The details of this effect as they appear in the non- 
relativistic SkI3 are shown quantitatively in Fig. U^ for 
selected neutron (left) and proton states (right), in both 
cases one level with large and one with small orbital an- 
gular momentum close to the Fermi energy. The upper 
panels show the radial density distributions Anr^ \4'{'i^)\ 
of the 2g and 3d neutron states and 2/ and 3p pro- 
ton states, where (j){r) is the radial component of the 
single-particle wavefunction ip{r). The radial density 
is shown for the state with larger total angular mo- 
mentum only. The middle panels shows the integrand 
rWg,r {i + ^m + lMe+i/2\^ -{i- kMe-i/2\!] which 
enters the calculation of the spin-orbit splitting (^) , while 
the radial component of the spin-orbit potential W^ is 
shown in the lower panels. Besides the familiar attrac- 
tive peak at the surface of the nucleus, the central de- 



pression of the density leads to a repulsive peak of the 
spin-orbit potential around r « 3 fm. The total spin- 
orbit splitting now depends sensitively on the location 
of the radial wavefunctions. The neutron 3d and proton 
3p states with three nodes but small angular momentum 
have large overlap with both the repulsive and the attrac- 
tive part of the spin-orbit potential (note that small radii 
are supressed only with 1/r and not as usual with l/r^), 
leading to nearly vanishing spin-orbit splitting, while the 
neutron 2g and proton 2/ states with only two nodes 
feel only the spin-orbit potential at the nuclear surface 
(and have much larger overlap with this than the small- 
angular-momentum states), showing the usual spin-orbit 
splitting. 

Note that this is a polarization effect that is natu- 
rally included in the self-consistent description of nuclei 
but cannot occur in semi-microscopic approaches like the 
"Extended Thomas- Fermi- Strutinski Integral" method 
(ETFSI) 1^,0 or macroscopic-microscopic models |6^] 
with prescribed densities and/or single-particle poten- 
tials, where one has a very restricted variational free- 
dom of the density profile only (ETFSI) or no degree-of- 
freedom in the density distribution and single-particle po- 
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FIG. 12. Distribution of the mass density from spherical calculations with Skl3 in the region of the Z = 120, A'^ = 172, and 
N = 184 shells. 



tentials at all (macroscopic-microscopic). Looking at the 
spectrum calculated with the FY model, the spin-orbit 
splitting is indeed much larger than in self-consistent 
models, especially for the 3p proton and 2g neutron states 
which are crucial for the A^ = 172 shell closure. Compar- 
ing Fig. g with Fig. |j for ^^^114 one immediatley sees that 
the change in the single-particle spectra of both protons 
and neutrons predicted by FY is much smaller when go- 
ing from 298x14 to ^^^120 than in all self-consistent mod- 
els. 

Figure O shows the profile of the total density in 
even-even nuclei in the region of the Z = 120, A^ — 172 
and N = 184 shells as they result from spherical calcu- 
lations with SkI3. This demonstrates that the density 
profiles are coupled to the shell closures (and vice versa). 
At large neutron numbers N > 184 all nuclei have the 
usual density profiles, while going below N = 184 the 
nuclei immeadiately show a central depression that is 
most pronounced for nuclei with Z = 120. It is note- 
worthy from Fig. n2^ that the central depression of the 
density distribution is coupled to the neutron number 
- it disappears for all neutron numbers above N = 184, 
while the density profiles of nuclei with constant neu- 
tron number but different proton number look very sim- 
ilar. The reason for this is that the last filled neutron 
levels below the N = 172 gap - 2gQ/2+ and IJ15/2- and 
2<?7/2+ ~ all have large orbital angular momentum and are 
therefore mainly located at the nuclear surface. Going 
from N = 172 to N = 184 only levels with small angular 
momentum - 8^5/2+, 8^3/2+ and 4si/2+ ^ are occupied 
which have a large probability distribution at small radii. 
This means that the unusual density distribution of nu- 
clei around 17I12O is simply caused by the filling of the 



neutron levels which have the same ordering in all mod- 
els investigated here. This effect thus should occur in 
non-self-consistent models as well. And indeed the den- 
sities calculated from the FY model (plotted in the upper 
panel of Fig. O) show the same behaviour as the densi- 
ties from the self-consistent models, although the effect 
is weaker here. But unlike the non-self-consistent models 
with prescribed potentials, the densities in self-consistent 
models are fed back into the potentials which amplifies 
the effect by driving the wavefunctions to larger radii. 
Additionally the self-consistent spin-orbit potentials are 
influenced which in turn causes the Z — 120 proton shell 
closure. 

The same effect which creates the Z = 120 proton shell 
is responsible for the appearance of a magic neutron num- 
ber N = 172. The gap at iV = 172 depends sensitively on 
the amplitude of the spin-orbit splitting of the id neutron 
levels above this gap. Therefore it occurs again only for 
the RMF parametrisations and the generalized Skyrme 
functional Skl3. It can be expected that this neutron 
shell closure is restricted to nuclei with a prominent cen- 
tral depression of the density like the Z = 120 proton 
shell closure. Fig. |3| shows the single-particle energies 
of the neutrons in the chain oi N = 172 isotones calcu- 
lated with Skis. The N = 172 gap is largest for Z = 120 
in agreement with our findings for the 62n in | pl| . Al- 
though all these A^ = 172 isotones show a central depres- 
sion of the density distribution, for those those around 
Z = 120 the decrease in density when going to small radii 
is steepest. This gives the largest (positive) peak in the 
spin-orbit potential and therefore the smallest spin-orbit 
splitting of the neutron 3d levels which in turn gives the 
largest gap in the spectrum. 
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FIG. 13. 1 Single-neutron levels in the vicinity of the Fermi 
energy for the isotones of A'^ = 172 versus the proton number, 
computed with SkI3. 



F. The Shell Structure of ?82l26 

The question weather Z — llA oi Z — 126 is the next 
spherical shell closure beyond the experimentally known 
Z = 82 is as old as the first extrapolations of nuclear shell 
structure to superheavy nuclei in simple models. While 
Z ^ 126 corresponds to the largest experimentally known 
magic neutron number, Z = llA has no counterpart for 
the neutrons. A large number of self-consistent models 
predict Z — 120 to be the next proton shell closure, but 
there are some parametrisations predicting Z = 126 as 
an alternative. 

Figure |lj shows the two-proton shell gap S2p for the 
chain of Z = 126 isotopes calculated with the forces as 
indicated. For SkP and SkM* two Skyrme forces forces 
which both have a large effective mass this is a major 
spherical shell closure. As in case oi Z — 120 the shell 
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FIG. 14. Two-proton gap in the chain oi Z = 126 isotones 
calculated with the parameterizations as indicated. 




FIG. 15. The same as in Fig. 



but for ?84l26. 



closure is neutron- number dependent, it fades away when 
going to neutron numbers beyond N = 184. For most 
other Skyrme forces there is only a slight enhancement 
of ^2p in a small vicinity around TV = 184 which cannot 
be interpreted as a shell closure. The forces with "rela- 
tivistic" spin-orbit coupling, i.e. all RMF forces and SkI3, 
predict very small shell gaps only. 

This is reflected in the single-particle spectra, see 
Fig. UB. Contrary to the appearance of the Z = \\A and 
Z = 120 shell closures, which can be explained simply 
by looking at the spin-orbit splitting of adjacent proton 
levels, the situation is more complicated for the Z = 126 
shell closure. The liii/2+ proton state which lies above 
the Z — 126 gap is widely separated from the deeply 
bound l«i3/2+ state. Therefore the appearance of the 
magic number Z = 126 does not depend only on the am- 
plitude of spin-orbit splitting but on the relative distance 
of levels with different orbital angular momentum as well, 
although all relativistic forces with overall small spin- 
orbit splitting show no shell closure at Z — 126. Remem- 
bering that states with large angular momentum have 
systematically too small single-particle energies and that 
the spin-orbit splitting predicted by the standard Skyrme 
forces and SkM is too large in heavy nuclei - both would 
reduce the Z — 126 gap - the occurrence of a proton shell 
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closure at Z — 126 is very questionable. 

Comparing the single-proton spectra of 1I4II4 (Fig. |j) 
and 1I4I26 one sees immedeatly that the gap aX Z ~ 126 
becomes much smaller with increasing proton number. 
An exception is the non-selfconsistent FY model, here 
the relative distances of all proton and neutron have only 
slightly changed. This gives a further example for the 
strong dependence of the shell structure of superheavy 
nuclei on the nucleon numbers in self-consistent models. 

For all forces the Fermi energy is positive which means 
that 184 126 is predicted to be unstable against proton 
emission. However, owing to the large Coulomb barrier 
in superheavy nuclei we expect that this nucleus decays 
trough other more common channels. 



G. Spin-orbit splitting in superheavy nuclei 

We have seen that the predictions of self-consistent 
models for the spin-orbit splitting in superheavy nuclei 
shows a pronounced dependence on the nucleon numbers 
and the orbital angular momentum of the single-particle 
states. This is summarized in Fig |l^. The upper panel 
shows the spin-orbit splitting of the 3/ (white markers) 
and 3p (black markers) proton states, while the lower 
panel shows the splitting of the 2g (white markers) and 
3d (black markers) neutron states in the nuclei as indi- 
cated for all forces under investigation. The trivial trend 
with the orbital angular momentum (. of the states is re- 
moved dividing by 2£ -|- 1 , see Eq. (^) . 

While in the non-self-consistent Folded- Yukawa (FY) 
model all states have nearly the same renormalized spin- 
orbit splitting, there are large differences between the 
self-consistent models. The predictions of the forces for 
certain states in certain nuclei differ as such, but there 
are clearly visible trends with nucleon number and or- 
bital angular momentum which occur for all parameteri- 
zations. Picking out one force, one sees in most cases the 
same pattern: The spin-orbit splitting of a given state 
in ^84126 is larger than in ^g4ll4, while it is smallest in 
17^120. The (renormalized) splitting of states with large 
orbital angular momentum is always larger than the split- 
ting of states with small orbital angular momentum. As 
already discussed above, this is related to the shape of 
the nuclear density distribution and the effect is largest 
in iy|l20, for which most self-consistent forces predict a 
pronounced central depression in the density. 

There is a difference between protons and neutrons. 
While the splitting of the 2g neutron state is compara- 
ble in all nuclei (although it follows the trend mentioned 
above), the differences with mass number for the 2/ pro- 
ton states is much more pronounced. 

There are large differences between the various forces. 
The parameterizations can be divided into three groups 
which differ in the isospin dependence of the spin-orbit in- 
teraction: standard Skyrme forces (SkP-Skll), extended 
Skyrme forces (SkI3, SkI4) and RMF forces (NL3, NL- 
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FIG. 16. Amplitude of the spin-orbit splitting of sev- 
eral superheavy nuclei as predicted by the mean-field forces 
as indicated. The Is-splitting is weighted with l/(2£ + 1) to 
remove the trivial dependency on the orbital angular momen- 
tum. This shows nicely that in self-consistent models the 
is-splitting has an additional state-dependency that does not 
occur in simple potential models like FY (In the modified os- 
cillator model the splitting is simply the k parameter in the 
potential) that is related to the shape of the density distribu- 
tion. 



Z, NL-Z2, and NL-VTl). The standard Skyrme forces 
in most cases predict larger spin-orbit splittings than 
the RMF forces. Like in the case of the known nuclei, 
the predictions of the extended Skyrme forces SkI3 and 
SkI4 do not stay in between the predictions of standard 
Skyrme forces and the RMF. Again, SkI3 gives much 
larger spin-orbit splittings than the RMF forces with a 
similar isospin dependence of the spin-orbit interaction, 
while SkI4 stays in betwen standard Skyrme forces and 
RMF forces for neutrons, but gives the largest splittings 
for proton levels. For SkP, the force with large effective 
mass m* /m = 1.0 and the smallest spin-orbit parameter 
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64, the results look somewhat different as it wa s alre ady 
the case for the known nuclei discussed in Sect. |VB| , the 
spin-orbit splitting of the large angular momentum states 
and the dependence of the amplitude of the renormalizcd 
spin-orbit splitting on the orbital angular momentum is 
smaller than in other standard Skyrme forces. 

The predictions for shell closures are sensitive on the 
isospin dependence of the spin-orbit interaction and the 
isoscalar effective mass. But there are additional depen- 
dencies of the spin-orbit splitting than the already men- 
tioned ones as can be seen by comparing SkI3 and the 
RMF forces, which have similar effective mass and isospin 
dependence of the spin-orbit interaction. 



VI. SUMMARY AND CONCLUSIONS 

We have investigated the influence of the isospin de- 
pendence of the spin-orbit force and the effective mass on 
the predictions for spherical shell closures in superheavy 
nuclei. 

We have introduced two new RMF forces: NL-Z2 and 
NL-VTl, both employing the standard nonlinear ansatz 
for the Lagrangian, but NL-VTl complemented with ten- 
sor couplings of the isoscalar and isovector vector fields. 
Both are fitted to the same set of experimental data 
as the recent Skyrme parameterizations Sklx. The ten- 
sor coupling changes the relative distances between the 
single-particle states, but it has no visible influence on 
spin-orbit splittings in heavy and superheavy nuclei. 

To test the predictive power of the models, we have 
compared the experimental and calculated single-particle 
spectra in ^°*Pb, the heaviest known spherical doubly- 
magic nucleus so far. Already in this nucleus used in 
the fit of all forces investigated here we see large differ- 
ences between calculations and experiment and among 
the forces. States with large angular momentum are 
shifted to too small single-particle energies and none of 
the self-consistent models gives the proper level ordering. 

The predictions for shell closures are found to be sen- 
sitive to the isospin dependence of the spin-orbit interac- 
tion and the isoscalar effective mass. The uncertainties 
of these quantities in the description of smaller nuclei 
amplify when going to large mass numbers, making pre- 
dictions for superheavy nuclei a demanding task. 

The occurrence of proton shell closures in self- 
consistent models depends strongly on the neutron num- 
ber (and vice versa), even when looking at spherical nu- 
clei only. This effect can be explained in terms of single- 
particle spectra as a coupling of the spin-orbit field to 
the profile of the density distribution (of protons and 
neutrons separately) which undergoes dramatic changes 
in superheavy nuclei. This is an effect of self-consistency, 
it cannot occur in models where the density distribution 
has only a restricted degree of freedom like the semi- 
microscopic ETFSI approach or has even no degree-of- 
freedom at all like in case of macroscopic-microscopic 



models. In the region around 172 120 all forces with 
small effective mass predict a deep central depression of 
the nuclear density, which induces an unusual shape of 
the spin-orbit potential that causes an additional state- 
dependence of the spin-orbit splitting. In some cases the 
usual level ordering of spin-orbit coupled states is even 
reverted. 

The change of the single-particle spectra of both pro- 
tons and neutrons when varying proton and neutron 
number is much larger in all self-consistent models than 
in non-self-consistent approaches, which was shown on 
the example of the Folded- Yukawa model. 

The only self-consistent force which predicts Z — llA 
for the next spherical magic proton number is the ex- 
tended Skyrme force SkI4. Although SkI4 gives a very 
good description of the binding energies in known (de- 
formed) superheavy nuclei ||32| and reproduces the kink 
in the isotopic shifts of the mean-square radii in heavy 
lead nuclei, it overestimates the spin-orbit splittings of 
proton states in heavy nuclei by 60 — 80%. This discrep- 
ancy between this very good description of bulk proper- 
ties and a rather poor description of details of the single- 
particle spectra is yet to be understood. Since a possible 
proton shell closure at Z = 114 is caused by a large spin- 
orbit splitting, the unique prediction of SkI4 is very ques- 
tionable. On the other hand, all RMF forces, which are 
in very good agreement with experimental data for spin- 
orbit splittings throughout the chart of nuclei predict a 
magic Z =120. 

In summary this gives a strong argument that the next 
magic proton number is Z = 120, coupled with a magic 
neutron number N = 172, still a far way to go from the 
heaviest presently known nucleus ^el US- 
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APPENDIX A: DETAILS OF THE MEAN-FIELD 
MODELS 

1. The Skyrme Energy Functional 

The Skyrme energy functionals are constructed to be 
effective interactions for nuclear mean- field calculations. 
For even-even nuclei, the Skyrme energy functional used 
in this paper 



£ == fkin[T] + i?sk[p, T, J] -f £c{pj\ ~ £c 



(Al) 



is composed of the functional of the kinetic energy fkin, 
the effective functional for the strong interaction £sk and 
the Coulomb interaction Eq including the exchange term 
in Slater approximation, and the correction for spurious 
center-of-mass motion £cm- The energy functionals are 
the spatial integrals of the corresponding Hamiltonian 
densities H. 

£[p,r,J]^ Jd\-H[p{r),T{v),3{r)] . (A2) 



The actual functionals are given by 



'Wkin - 


2m 




Hc-- 


_e^ /• 3 Pp(r)pp(r') 36^/3^ 
2 J |r-r'| 4 Vtt^ 


Hsk = 


-'ip' + b^pr-'-lpAp+'fp'^^^ 




-E 


2 T / 2 A 
JPq+ hPqTq " 1^ Pq^Pc 




+Wls 





1/3 



Pp- 



4/3 



(A3) 

(A4) 



"3 Q 2 

^P Pq 



(A5) 



with various possibilities for the spin-orbit interaction 



n 



;std) 

LS 



''•LS ^ ''•LS 



-T/(ext) 



bi(p-^ -S + ^PqV ■:Sq 

q 

q 

^4 PV • J - 64 2J Pq V • 3q 



(A6) 

(A7) 
(A8) 



'^LS ^^ reproduced from Hl™ setting b'^ = b4. 

The local density Pq, kinetic density Tq and spin-orbit 
current 3q entering the functional are given by 



Pq= J2 ^k IV'feP , 

Tq^ Y,vl\V^u? , (A9) 

ken, 

J, = -f ^ Ufe Vi V X (T-0fc - (V X CTV'fc)'^V'fc 



with q £ {p,n}. Densities without index denote to- 
tal densities, e.g. p — pp + pn- The -0^ are the single- 
particle wavefunctions and v'^ the occupation probabili- 
ties calculated taking th e res idual pairing interaction into 
account, see Appendix A3. The parameters bt and b'^ 
used in the above definition are chosen to give a most 
compact formulation of the energy functional, the cor- 
responding mean-field Hamiltonian and residual interac- 
tion p8|. They are related to the more commonly used 
Skyrme force parameters ti and Xi by 

5o = to(l + ixo) , 

b'l ^ iM^ + ^1) - t2{^ + ^2)] , 
^2 = i [3^1 (1 + 52:1) -f2 (1 + 5X2)] , 

b'2^^[Sh{^+Xi)+t2{^+X2)] , 

63 = ii3(l + ia;3) , 

b's^-Mh + xs) , 

c\ = -|(tia;i +t2a;2) , 

c; = -i(ti-f2) , (AlO) 

The actual parameters for the parameterizations used in 
this paper are summarized in Table O. 

The single-particle Hamiltonian is derived variation- 
ally from the energy functional. One obtains 



(All) 



hq = -\I ■ BqV + Uq - 


- iWq • V X (7 


with the mean fields 




^q ^ c ' '-'q = X 

dTq dpq 


w - ^^ 

' "-SJq 



(A12) 



For all forces, a center-of-mass correction is employed. 
For the Sklx and SLyx forces it is calculated perturba- 
tively by subtracting 



£c 



2mA^ ™^ 



(A13) 



from the Skyrme functional after the convergence of the 
Hartree-Fock equations, while for SkM* and SkP only 



kefiq 



the diagonal direct terms in (A13) are considered self- 
consistently in the variational equation |4^. For all but 
the SLyx forces this is the procedure used in the or igina l 
fit. For SLy6 and SLy7 the microscopic correction (A13) 
was considered in the variational equations and therefore 
giving a contribution to the single-particle energy. How- 
ever, for large nuclei as discussed here the contribution of 
( A13| ) to the single-particle energies is negligible because 
the matrix elements are weighted with l/(2mA) com- 
pared to the contributions from the energy functional. 
We have therefore omitted this feature and follow the 
suggestion of [ ^ to use the perturbatively calculated 
correction from (A12) instead. 
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TABLE II. Parameters of the Skyrme energy interactions. The ti, Xi, 64, 6i and a are the parameters of the Skyrme 
functional (AS), h^ /2m is the constant in the calculation of the kinetic energy (A3). 



Parameter 
to [MeV fm^] 
ii [MeV fm^] 
£2 [MeV fm^] 
t3 [MeV fm^+"] 

Xo 
Xl 
X2 
X3 

64 [MeV fm^] 

&4 [MeV fm^] 

a 

h^/2m [MeV fm^ 



SkM* 



SkP 



Skll 



SkI3 



SkM 



SLy6 



SLy7 



-2645.0 


-2931.70 


-1913.619 


-1762.88 


-1855.827 


-2479.50 


-2480.80 


410.0 


320.662 


439.809 


561.608 


473.829 


462.180 


461.290 


-135.0 


-337.41 


2697.594 


-227.090 


1006.855 


-448.610 


-433.930 


15595.0 


18708.97 


10592.267 


8106.2 


9703.607 


13673.0 


13669.0 


0.09 


0.29215 


-0.954536 


0.3083 


0.405082 


0.825 


0.848 


0.0 


0.65318 


-5.782388 


-1.1722 


-2.889148 


-0.465 


-0.492 


0.0 


-0.53732 


1.287379 


-1.0907 


1.325150 


-1.0 


-1.0 


0.0 


0.18103 


-1.561421 


1.2926 


1.145203 


1.355 


1.393 


65.0 


50.0 


62.130 


94.254 


183.097 


61.0 


62.5 


65.0 


50.0 


62.130 


0.0 


-180.351 


61.0 


62.5 


1/6 


1/6 


0.25 


0.25 


0.25 


1/6 


1/6 


20.733983 


0.733983 


20.7525 


20.7525 


20.7525 


20.73552985 


20.73552985 



TABLE III. Parameters of the RMF forces used in this investigation. The mass of the isovector vector-field rUp 
is not fitted and is the same for all forces. 



763 MeV 



Force 



(MeV) 



m^ 
(MeV) 



m^ 
(MeV) 



9p 
(fm-i) 



62 



63 



3 


939.0 


508.194 


782.501 


10.2170 


12.8680 


4.47400 


-10.4310 


-28.8850 


-Z 


938.9 


488.67 


780.0 


10.0553 


12.9086 


4.84944 


-13.5072 


-40.2243 


-Z2 


938.9 


493.150 


780.0 


10.1369 


12.9084 


4.55627 


-13.7561 


-41.4013 


-VTl 


938.9 


484.307 


780.0 


9.81307 


12.6504 


4.63432 


-13.2808 


-38.0773 



fu^/gu. = -0.102703 
fjgp = -4.71143 



2. Relativistic Mean-Field Model 

For the sake of a covariant notation, it is better to 
provide the basic functional in the relativistic mean-field 
model as an effective Lagrangian C. For the present ver- 
sion of the RMF used in this study, we can summarize it 
as 



C 



RMF 



^U 



-M 



-NM 



-'nonl 



(A14) 



where £n is the free Dirac Lagrangian for the nucleons 
with nuclcon mass ton , equally for protons and neutrons 



C-^ = ^vl ipk {i^fid^^ - "t-n) ipk 



(A15) 



fcesi 



The Lagrangians of the fields and their couplings to the 
nucleons are given by 



i(9^$,,, - a.$,,^) • a^$; - to5$^,^ • $ 









(A16) 



The model includes couplings of the scalar-isoscalar 
(<i>cr), vector-isoscalar {^ui,fj.), vector-isovector (^p,^), 



and electro-magnetic (A^) field to the corresponding 
scalar-isoscalar (ps), vector-isoscalar (p'^) and vector- 
isovector {p'^) densities of the nucleons as well as the 
proton density p^, which are defined as 



Y^ vl iJkipk, 
ken 

Y vl ■4>kip'4>k 



ken 



= Y^l ^pkfjf,^pk, 



ken 



E 



Pp = y . ^k Aip^k- 



(A17) 



lAaa is the nonlinear selfinteraction of the scalar-isoscalar 
field. All forces used in this paper employ the standard 
ansatz |3q| 









(A18) 



In case of the parameterset NL-VTl also a tensor cou- 
pling between the nucleons and the vector fields is con- 
sidered, which can be written as 



u 



-NM 



2ton 
with the densities 



*-,mP" + 



Z' -J- ^P $^ 



2ton 



■PC 



(A19) 
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p'i^d^Y. ^l^kT^^^k, 



ken 



Pt 



= d,Y.vl^Pk<j''''T^Pk, 



(A20) 



ken 



where a^'^ — {i/2)[y^,'j'^]. The masses rrii and coupHng 
constants of the fields are the free parameters of the RMF 
which have to be adjusted to experimental data. The 
actual parameters of the parameterizations used here are 
given in Table [II . The equation of motion of the single- 
particle states is derived from a variational principle 

efcToV-fc = [ - «7 • V + mN + 5 + 7^vH V-fc (A21) 



where S = g^^a and V^ = 3^*^.^ + ^P^P.fj. ■ ■? + 
^eA^{l + To) are the scalar and vector field respectively. 
A more detailed description of the model can be found 

in d]. 

For the residual pairing interaction and the center-of- 
mass correction the same non-relativistic approximation 
is used as in the SHE model, for NL-Z, NL-Z2 and NL- 
VTl by subtracting pcrturbatively the full microscopic 
correction (A13), while for NL3 the harmonic oscillator 



estimate Ec 



4:1A ^/^ MeV is subtracted as done 



in the fit of these parameter sets. 



3. Pairing Energy Functional 

Pairing is treated in BCS approximation using a delta 
pairing force |73,p0], leading to the pairing energy func- 
tional 



p . — _ 

'-'pair — , 



E 

= {P.n} 



v„ 



rXn 



(A22) 



where Xq = -'^J^ken^yo hukVk \'4'k? is the pairing den- 
sity including state-dependent cut-off factors /^ to re- 
strict the pairing interaction to the vicinity of the Fermi 

TABLE IV. Pairing strength Vn for the neutrons and Vp 
for the protons for the mean-field forces used in this study. 
m* /m is the isoscalar effective mass in infinite nuclear matter. 
Note that the absolute value of the pairing strength decreases 
with increasing effective mass. 



Force 


m* jm 


K [MeV fm^] 


Fp [MeV fm3] 


SkM* 


0.789 


-276 


-292 


SkP 


1.0 


-241 


-265 


Skll 


0.693 


-320 


-305 


SkI3 


0.574 


-340 


-351 


SkI4 


0.650 


-310 


-324 


SLy6 


0.689 


-308 


-320 


SLy7 


0.688 


-308 


-320 


NL3 


0.595 


-329 


-342 


NL-Z 


0.583 


-349 


-351 


NL-Z2 


0.583 


-343 


-350 


NL-VTl 


0.600 


-340 


-346 



surface [Q . u^ is the occupation probability of the given 
single-particle state and u^ = 1— u^- The strengths V-p for 
protons and T4 for neutrons depend on the actual mean- 
field parameterization. They are optimized by fitting for 
each parameterization separately the pairing gaps in iso- 
topic and isotonic chains of semi-magic nuclei throughout 
the chart of nuclei. The actual values can be found in Ta- 
ble 1^. The pairing-active space Q.q is chosen to embrace 
one additional shell of oscillator states above the Fermi 
energy with a smooth cutoff weight, see Mj for details. 



4. The Folded- Yukawa Single-Particle Potential 

We present here only the details needed for our dis- 
cussion. A more detailed discussion of the parameter- 
ization of the potentials can be found in |69[| and ref- 
erences therein. The single-particle Hamiltonian of the 
Folded- Yukawa single-particle model has the same struc- 
ture as the one of the Skyrme-Hartree-Fock model (All), 
but instead of calculating the potentials self-consistently 
from the actual density distributions, a parameterized 
guess for the functional form of the potentials is used. 
The nucleons have an effective mass of m* /m = 1 with- 
out any radial dependence, therefore B is simply given 
hy B = h /2m. The single-particle potential U is cal- 
culated from the folding of a Yukawa function with the 
sharp nuclear surface 



Uq{v) 



-J^ d\'— - 

Atto^ Jy |r — r'l/a 



(A23) 



where the integration is performed oder the nuclear vol- 
ume. Finally, the spin-orbit potential is given by the 
derivative of the nuclear potential 



(A24) 



W,{v) = -\{A)[ — ] VU, 



with the coupling constants Ap = 28.0 + 6.0 A/24Q and 
Ap = 31.5 -F 4.5 A/240. 
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